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I. INTRODUCTION 



The emergence of manganites as a technologically important material due primarily to the 
occurrence of colossal magnetoresistance (CMR) in e.g. Lai_ x Ca x Mn0 3 Jl| has rekindled 
the interest of the condensed matter physics community in the double-exchange mechanism. 
The physics of double-exchange has successfully correlated ferromagnetism and metallicity 
in the doping range of 0.2 < x < 0.4. The basic ingredient of double-exchange, proposed 
by Zener fifty years ago ||, is encapsulated within a simple two site model of Anderson 
and Hasegawa ^ [|, The latter has led to a plethora of theoretical and experimental 
investigations in recent years || [Tj, which have gone on to add one or the other feature to 
the original model, often without abundant care, as we shall argue here. 

The undoped manganite system is characterised by the presence of an incomplete d-shell 
of the Manganese ions, which consists of 3 electrons in the t 2g state and one electron in the 
e g state. The t 2g spins are deep inside the d-level and are assumed to be unimportant in the 
process of charge transfer. However, these spins do show a tendency to align antiferromag- 
netically both in the parent compound and in the completely doped system. The Coulomb 
energy cost for the e g electrons to hop onto the adjacent Manganese ion in the absence of a 
hole is very large. The undoped system is thus an insulator. 

We will henceforth refer to the e g electron as the itinerant electron, the t2 g electrons as 
the core electrons and the total spin in the t2 g level as the core spin. Doping the system 
with a divalent atom like Ca or Sr leads to creation of holes in the e g level in a fraction of 
the Manganese ions. The itinerant electron from the occupied e g level of one Manganese site 
can then hop into its nearest neighbor Oxygen site, which facilitates hopping to the nearest 
neighbor unoccupied e g level of another Manganese ion, thus leading to a finite conductivity. 
This process, called the "double exchange mechanism", plus the presence of strong Hund's 
rule coupling between the core spin and the itinerant electron, results in an indirect coupling 
between neighboring core spins . This in turn relates the magnetic order of the underlying 
lattice with the kinetic energy of the itinerant electron. 
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The conclusions derived by Anderson and Hasegawa can be summarised as follows. In 
the limit when Hund's coupling is infinitely strong, the itinerant electron would like to have 
its spin aligned with the local core spin. Additionally, if the core spin is treated classically, 
the appropriate axis of quantization is the direction of the core spin vector, as far as the 
itinerant electron is concerned. Now, as the latter hops, it has to readjust its spin to be 
realigned with the new core spin partner, amounting to a rotation of the quantization axis 
by an angle 9, which is the polar angle between the core spins Si and S 2 - From the property 
of the spin-1/2 rotation operator it follows that the hopping or the overlap matrix element 
t will be renormalized to tcos((9/2), assuming azimuthal symmetry. It then follows that 
if 9 equals tt, the core spins have antiferromagnetic (AFM) coupling and hopping of the 
itinerant electron is totally inhibited. On the other hand, if 9 equals zero, the core spins 
have ferromagnetic (FM) coupling and hopping is accentuated, thus synergising transport 
with ferromagnetism. This is the first result of Anderson and Hasegawa. The latter then 
proceeded to a quantum treatment of the core spins, but still operating within the infinite 
Hund's coupling limit. Interestingly, it turns out that the energy eigenvalues are identical 
to the earlier classical case, provided cos(6 l /2) is identified as a parameter which equals 
(Sq + l/2)/(2S + 1), where S = \ Si \ = \ S 2 | and So = | Si + S 2 + |, o being the spin 
of the itinerant electron (| a \ = 1/2). Curiously, since Si, S 2 and a add up to (25* + 1/2) 
in the FM case for large Hund's coupling, the parameter cos(#/2) would indeed be equal to 
unity, as in the classical case. But in the AFM case, the parameter reduces to 1/(2S + 1), 
which goes to the classical value of zero only when S —>■ 00, thus necessitating an additional 
constraint on the core spin, if the "classical" interpretation is to be taken seriously. Needless 
to say, in between FM and AFM cases, the parameter cos(#/2) would go through a set of 
discrete (and not continuous) values, thus pointing to the need of a more careful treatment 
of the core spins in the quantum case. 

One of the directions in which the Anderson- Hasegawa treatment has been extended is 
to recognize the importance of an additional superexchange term between the core spins 
proportional to Si.S 2 . It has been assumed in the literature till now that the superexchange 
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term can be simply taken as an extra term to be added to the Hamiltonian and that the 
large Hund's rule coupling affects only the process of charge transfer in these systems. In 
this paper we show inter alia that superexchange is itself modified in a nontrivial manner 
if the core spins are dealt with quantum mechanically, which leads to a correlated diagonal 
disorder in these systems, even in the cleanest samples. Our starting point therefore is 
the Hamiltonian for a two site one electron model, including the superexchange interaction, 
given by 

2 

H = -tJ2(4 T c 2r + h.c) - Jn^SiA + JSi.Sz. (1) 

T 1=1 

Here t is the hopping matrix element for the itinerant electron between the two sites, c\ T 
(c iT ) is the creation(annihilation) operator of the itinerant electron at site i having spin 

— * 

projection r, J H is the Hund's rule coupling strength, Si is the core spin at site i and <7; 
is the spin of the itinerant electron at the site %. The parameter J is the superexchange 
interaction strength between the core spins in the nearest neighbor sites. For our case we 

— * 

consider | Si \— S, i.e. the core spins on all the sites are taken to have the same value. 

With the preceding background, the motivation behind our work and the plan of this 
paper are as follows. We reiterate that the two site, single electron, double-exchange model 
is the simplest basic framework for interpreting a large number of fascinating properties of 
manganites. While Anderson and Hasegawa did provide a quantum solution to the model, 
especially in the limit of large Hund's rule coupling, subsequent authors seem to have gone 
ahead in a somewhat cavalier fashion, in our opinion, about the classical limit of infinitely 
large core spins. This has led to some confusion about the interpretation of parameters in 
the model which needs to be cleared. As the value of the core spins in most studied CMR 
systems is indeed finite — three-halves for manganese — it is important to delineate the 
quantum versus classical effects, especially while considering additional phenomena, e.g., 
polaron-induced hopping and thermodynamics. With this aim in mind, we organize the 
paper as follows. We present in Sec. II, the exact quantum mechanical solution to the 
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Anderson-Hasegawa model with an additional superexchange term, take the large Hund's 
rule coupling limit and discuss the outcome of a hitherto ignored 'site-disorder' term. In 
Sec. Ill we reexamine the issue of polaron-assisted hopping in the light of our fully quantum 
calculation. The results in this section are then employed in Sec. IV for the computation 
of heat capacity and phase diagram, wherein we also specify the difference between our 
results and those which treat the core spins classically. Finally, in Sec. V, we present some 
concluding remarks. 

II. EXACT SOLUTION FOR THE ANDERSON-HASEGAWA MODEL 



In order to find the ground state of the system we follow the quantum mechanical 
calculation carried out by Anderson and Hasegawa ||. We first note that the Hund's 
rule coupling term, proportional to Jh is diagonal in the states given by | ipf) = 
| Si, ± \),S 2 ;S Q ,M > and | ipf) = I Si,~,S 2 ,(S 2 ± \);S G ,M >, while the hop- 

ping part of the Hamiltonian connects these two sets of states, corresponding as it does to a 
recoupling of the itinerant electron's spin (1/2) from the site spin Si to 5*2 and is thus given 
by the Wigner 6j (or Racah) coefficient (W)|§. Here, M = Sq. The superexchange term 
proportional to J, is off-diagonal in the basis states chosen above. However, it is diagonal 
in the states given by | 0(5")) = | ^, Si, S 2 , (5"); S , M), where S' = \ Si + S 2 |. We can then 
relate the states | 4>f 2 ) to the states | 0(5")) through appropriate Racah coefficients again 
and we find, 




#+1) W (lSiS'S 2 ; (Si + l)S') | 0(5')) 



vW+1) W (^SiS'S 2 ; (5 2 + 1)5^ | 0(5')). (2) 



Clearly, since 5" (the total core spin) must couple to the itinerant electron spin to give the 
total angular momentum 5"o, the only values of 5" to be summed over are 5" = 5o + | and 
5o — 2 ■ The particular Racah coefficients which occur (with Si = S 2 = S), have convenient 
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closed expressions such as, 



= cos(a/2) | <P(S - i)> + sin(a/2) | 0(S O + ±)>, 
^f) = -sin(a/2)|0( 1 So-i)) + cos(a/2)|0( 1 S o + ^)), 



(3) 



where, 



cos(a/2) 



(25 + S + |) 



1/2 



(4) 



2(25+1) 

The relations between | ijjf) and | (f)(S ± |)) states are the same as those between | ipf) 
and | 4>(Sq ± i)). Note that the expression for cos(o;/2) can be written in terms of cos(6*/2), 
which actually gives the relation between a and 6 as a = 9/2. Thus the Hamiltonian matrix 
in the space of | ipf) and | ipf) can be written as follows : 

/ Pi Pi -tcos(0/2) -tsin(0/2) N 



where, 



*sin(0/2) -tcos(6/2) 



Pi 



-tcos(0/2) tsin(0/2) 



\ -tsin(0/2) -tcos(9/2) 



P 1 = 3 [R x sin 2 (a/2) + i? 2 cos 2 (a/2)) - —5, 

P 2 = -(i?! - J R 2 )cos(a/2)sin(a/2), 

P 3 = ^(i2iCOS 2 (a/2) + i2 2 sin 2 (a/2)) + ^(S' + l), 



(5) 



(6) 



and, 



Ri=[S + ^j (s + l) -2S(S + 1), 
R2=[S -l) (s + l)-2S(S+l). 



(7) 
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The eigenvalues (E) of the Hamiltonian matrix are obtained from, 



2E = ^ + K^J) ± ^At 2 + K 2 {J) + K 2 (J) ± Atcos(9/2)^K 2 (J) + K 2 (J), (8) 



where, 



K X {J) 

K T {J) 

K*{J) 
K 3 (J) 



J 



S + -) -2S(S+1) 



M2S+1) 



+ K 2 (J) 



J cos(a)(S'o + -) 
J sin(a)(5'o + -). 



(9) 
(10) 
(11) 



While the exact result of Eq.(^) may be of interest in its own right, we examine the limit 
of large Hund's rule coupling, by expanding the square root term (upto 0(1/ Jr))- We find 
that the lowest energy eigenvalues are given by 

E m = -M - icos(0/2) + - [S'(S' + 1) - 2S(S + 1)] + Jsin 2 (a/2)(So + i), (12) 

where S' = S - 1/2. 

The first two terms in Eq. (|12"D are the terms obtained by Anderson and Hasegawa. 
We emphasise once again that the parameters cos(0/2) and cos(a/2) take discrete values 
which depend on the quantum values of the core spins. The third term is the result of 
the superexchange interaction in the absence of the itinerant electron. The fourth term, a 
novel one, is purely due to the double-exchange mechanism in the presence of the itinerant 
electron. The explicit form of this term is given by, 

J2S-S'-, 



AE.i 



2 25 + 1 



•(S' + l). 



(13) 



It is to be noted that in the ferromagnetic limit (i.e. S' = 2S) AEj vanishes exactly. 
There are primarily three important points to be made about AEj : 
(a) This term is an on-site term : It does not involve physical transfer of the itinerant 

electron from one site to the other. 
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(b) This term vanishes in the absence of the itinerant electron. Thus, it exists only on 
the site at which the electron resides and hence, at any site i, it will be proportional to n^, 
where n« is the number of itinerant electrons at the site % and is taken to be either 1 or in 
the absence of double occupancy. 

(c) At any site i this term depends on 3 spin values : the spin of the itinerant electron 
on the site i (<7j), the core spin at site i (Si) and the core spin on the neighboring site (Sj), 
j being the neighbor of the site i. 

Thus, the extra energy term AEj corresponds to a site energy term which is correlated 
with the core spins of the nearest neighbors. We may therefore propose an effective double- 
exchange Hamiltonian in the full lattice as 

#eff = Yl e * n * - 1 E cos(^/2)(ct Cj + H.C.) + JijSi-Sj, (14) 

i <ij> <ij> 

where 

3 

and 

S[ 3 =| Si + Sj I . (16) 

We have explicitly taken in order to accomodate effects of anisotropic superexchange 
also. 

The classical limit of the extra term is given by 

AE? = ^-{2S + 1) [1 - cos(0/2)] cos(0/2), (17) 

which goes to zero in both the ferromagnetic as well as the antiferromagnetic limits. Again, 
we see a clear distinction between the classical and the quantum results in the antiferromag- 
netic limit. We emphasise that in taking the purely classical expression, one actually loses 
the effect of the quantum fluctuations which are present in these systems not only because 
of the fluctuating spins but also due to the on site disorder and the hopping, correlated 
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with the spins on the lattice. A variety of interesting physical phenomena could be studied 
by taking into consideration the quantum Hamiltonian. One of these concerns the polaron 
effects, which are discussed below in Section III. 



III. ANDERSON-HASEGAWA-HOLSTEIN MODEL 

Experiments in manganites - both thermodynamic and transport - seem to suggest the 
importance of polaron formation and the consequent localization of charge carriers ||. The 
minimal model which reflects such lattice carrier interaction on the double-exchange can be 
introduced by dovetailing the Holstein mechanism on the Anderson-Hasegawa Hamiltonian. 
Therefore, in view of the results presented in Sec. II, in the limit of large Hund's rule 
coupling, we may write a two site Anderson-Hasegawa-Holstein Hamiltonian as, 

2 g + 1 2 

H = e Y n ™ ~ Y t 2 °s + 1 ( c i° C2a + c 2^) + 9iUo Y Y ni °( bi + 6 *) 

i=l a a i=l a 

2 

+ g 2 oj [niAh + h\) + n 2a {b x + &{)] + lu ^ h\b, + JS X .S 2 + AE J} (18) 

a i=l 

where, g\{g%) denotes the on-site (intersite) electron-phonon coupling strength and e is the 
bare site energy. This site energy is to be contrasted with the site energy which was derived 
in the previous section. Here, since we are interested in doing a two-site model, the extra 
term in the Hamilonian due to the superexchange interaction, has been introduced simply 
as AEj. Note that we have considered a single phonon mode for interatomic vibrations of 
frequency uq for which 6, and b\ are the annihilation and creation operators. 

We separate out the in-phase mode and the out-of-phase mode by introducing new phonon 
operators a = (61 + 62)/ v^2 and d — (bi — b 2 )/ \pl in the Hamiltonian. The in-phase mode 
does not couple to the electronic degrees of freedom whereas the out-of-phase mode does, 

leading to a Hamiltonian Hd, given by, 

2 / S + - \ — - 

H d = u Q cfld + eJ2 n i-t ( 2 g + 1 ) ^ l ° 2 + h - c ^ + d-M^i ~ n 2 )(d + d ] ) + JSi.S 2 + AEj, 

(19) 



which represents an effective electron-phonon system. Following |10| we use a Modified 
Lang-Firsov (MLF) transformation with variable phonon basis and obtain, 



R tt „-R 



H d = e u H d e 

S + - 

= u dd + V e p ni - t-^ — \ [c\c 2 exp(2A(d f - d)) + c\d exp(-2A(d f - d))] 

+ uj (g^-X)(n 1 -n 2 )(d + d^ + jJ2Si-S J +AE Jj (20) 

<ij> 

where R = A(ni — n 2 )(<f~ — d), A is a variational parameter related to the displacement of the 
d oscillator, g_ = (gi — g 2 )/\/2 and t p = e — cu (2g_ — A)A. The basis set is given by |±, N) = 
-^(c\±c\) \0) e \N), where |+) and |— ) are the bonding and the antibonding electronic states 
and \N) denotes the A^th excited oscillator state within the MLF phonon basis. The diagonal 
part of the Hamiltonian in the chosen basis is treated as the unperturbed Hamiltonian 
(Hq) and the remaining part of the Hamiltonian H% = Hd — Ho, as the perturbation. The 
unperturbed energy of the state |±, N) is given by 



E^ N = (N,±\H \±,N) 

r \t 

(2A) 



Nu Q + e p =F t eff 



i=o l - 



+ JS 1 .S 2 + AE J (21) 



where t e ff = t 2S+1 ex P ( — 2A 2 ), Nd = rjyz^ ^ ■ The general off-diagonal matrix elements 
of Hi between the two states |±, N) and |±, M) may be calculated for (N — M) > as in 
Ref. |0|. 



The unperturbed ground state is the |+)|0) state and the unperturbed energy, = 
tp — t e ff + JSi.S 2 + AEj. However, in this exact quantum limit of core spins, for given 
values of g_ and J, E^ can have four values corresponding to ferromagnetic (FM), canted 
1 (CAl), canted 2 (CA2) and antiferromagnetic (AFM) orientation of the two spins for 
I 5*12 | = | 5*1 + 5*2 |= 3, 2, 1, respectively. Minimizing the unperturbed ground state energy 
A is calculated and is given by 
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We have evaluated the perturbation correction to the energy upto the sixth order and 
the wave function upto the fifth order. The convergence of the perturbation series is very 
good for t/u)Q < 1. To obtain the ground state spin order of the core spins we calculate the 
energy for each set of values of g_ and J with four possible Si% and find out the combination 
for which the energy is the minimum. Further, to study the effect of an external magnetic 
field (h) we include a term —cjhb(S\ + S-z).h to the Hamiltonian in equation fll8|), g being 



the Lande g factor. We assume that the external magnetic field is along the direction of S\2 
and is expressed in units of /x e //(= PA*b)=1- 

It is expected that the charge transfer from site '1' to '2' depends on the spin order of 
the core spins as well as the electron-phonon interaction. In the double-exchange model, the 
effective hopping reaches its maximum value in the ferromagnetic state and decreases as it 
approaches the antiferromagnetic limit. Moreover, in a lattice, the electron produces lattice 
deformations and which in turn localize the electron for strong electron-phonon coupling. 
To study the polaronic character one calculates the static correlation functions (n 1 -u 1 ) and 
(711112)0, where u\ and U2 are the lattice deformations at sites 1 and 2 respectively, produced 



by an electron at site 1 |10|, |TT|| . In the present report with a two-site one electron model, 



following [p| , we calculate — < 7ii(«i — u 2 ) >o / 9- = — g — and study the nature of the 



polaron crossover for different ranges of <?_ and J. In the 'large' polaron limit this parameter 
takes a small value, while with increasing electron-phonon coupling it tends to unity, showing 
a distinct crossover from 'large' to 'small' polaron behavior. The measure of derealization 
of the electron for various ranges of g_ as well as J will be evident from the kinetic energy. 
So we have also calculated the kinetic energy, given by, 



S + 1 

= -E Kin =< ^ G |t exp(2A(d t - d)) + eld exp(-2A(rf t - d)) }\tfj G >, 



(23) 
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where ipc is the ground state wave-function, is evaluated upto the fifth order in the pertur- 
bation. The numerical evaluation of E Kin will be presented below in Sec. IV. 



IV. PHASE DIAGRAMS AND SPECIFIC HEAT 



Recently, there have been many experimental reports on manganites at low doping and 



low temperatures with and without an external magnetic field [|12| , |13| , pHfl . Okuda et al 
have estimated the electronic specific heat for La\- x Sr : x MnOs in the ferromagnetic regime 
and concluded that the carrier mass-renormalization near the metal-insulator transition at 
x = 0.16 is minimal. They have also observed a decrease in the low temperature specific 
heat in the presence of a magnetic field. Motivated by these observations, we have carried 
out a calculation of the specific heat, based on the partition function of the system which, 
from a cumulant expansion upto the 2nd order, is given by ]T3|, 



Z{f3) = Z ((3)exp(- f V f dp H {HW)Hi{p"))), (24) 
Jo Jo 

where Z ((3) = Tr(e~ pHo ) ; #i(/3) = e^^e"^ , and (3 = The expression () denotes 

the usual canonical averaging. The specific heat is then calculated (in arbitrary units) from 
the well known relation: 



c " = '^ip lnzm ' (25) 

and in the low temperature regime, to which only the zero-and one-phonon states contribute. 

If the localized core spin at each site is | then the possible values of | Si + 02 \— S12 are 
3, 2, 1 and 0. The ferromagnetic (FM) and antiferromagnetic (AFM) orders are obviously 
related to 5*12 = 3 and 0, whereas S12 = 2 and 1 are referred as canted 1 (CAl) and canted 
2 (CA2) states respectively. The Fig. 1 shows the phase diagram for the four possible spin 
orders for our system, in the g_ vs J plane. For small values of g_ and J, the FM state is the 
most stable one, and with increasing J, the ground state first becomes CAl and then CA2. 
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For a very large value of the superexchange interaction J, the system is in an AFM order 
for any value of g_. However, with increasing electron-phonon interaction g_, the CAl and 
CA2 phase become narrower. Indeed, for larger values of g_ the FM state appears for very 
low J but with a small increase of J the system transits to the AFM phase. The CAl and 
CA2 phases in fact do not appear at all as the phase changes from the FM to AFM state 
with increasing superexchange interaction J, for large values of g_. It can be further shown 
that for a very large value of g_ the ground state is AFM for any value of J. 

It is evident from the phase diagram that for a particular value of J the ground state 
changes as the electron-phonon coupling (g_) increases (Fig. 1). But the change of phase 
from one to another is not continuous with g_, for the quantum consideration of the core 
spins. This is shown in Fig. 2. For small values of J = 0.01, the FM state exists even for 
a large value of g_ and then it sharply changes to the AFM state. On the other hand, for 
larger values of J = 0.04, 0.09, the system passes sharply to the canted phases (CAl and 
CA2) and then the AFM state, with increasing g_. For J =0.04 the CAl and CA2 regimes 
are very narrow and for J =0.09 the CAl and CA2 orders persist for a wider range of 
This is to be contrasted with the classical core spin model in which a similar study shows 
that the transitions to different core spin orientations are continuous for the same range of 
values for g_ and J []lT |. In the classical case only three phases(FM, AFM and Canted) are 
present. The relative angle 9 between classical core spins can take any value from to tt, so 
any spin orientation other than FM (9 = 0) and AFM [9 = tt) yields a canted phase. Hence, 
in the classical limit of the core spins, for certain values of J, the FM-AFM transition is a 
smooth and continuous transition with g_, whereas for spin |, the FM-AFM transition with 
g_ is never continuous for any J. 

The probability of hopping of the itinerant electron from site to site is a maximum in 
the FM state as would be expected from the double-exchange mechanism. But, for a very 
strong electron-phonon coupling g_, the electron may be localized forming a small polaron. 
For low values of J we find both small and large polaron ground states in the FM phase 
(Fig. 3). The large to small polaron crossover is indicated by the relative deformation of 
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the two lattice sites which is measured by the static correlation function In Fig. 3 

the kinetic energy is large for small values of where the polaron is large, and for 
large g_, the kinetic energy reduces rapidly, while rises, showing a smooth crossover to 
the small polaron regime. The classical and quantum formulations of the double-exchange 
model turn out to be the same in the FM limit of the core spins. So the nature of the kinetic 
energy and the polaron crossover in the FM state, as shown in Fig. 3, will be unaltered in 
the classical limit of core spins. However, in the AFM limit the two approaches (classical 
and quantum) are not equivalent, as has been argued earlier also. In the S — > oo limit the 
hopping probability is zero for the AFM case, while for S — | the parameter modifying the 
hopping probability t, takes a finite value 0.25 resulting in a finite charge transfer, even in 
the AFM limit. 

In Fig. 4 we show the nature of variation of the kinetic energy as well as the polaron 
crossover in different magnetic ground states. For J = 0.09 the ground state is FM for low 
<7_, and with increasing g_, the ground state changes sharply to CAl, CA2 and lastly to 
the AFM state. Since at each transition (from FM -> CAl -> CA2 -> AFM) the effective 
hopping reduces due to the double-exchange interaction, it is obvious that the kinetic energy 
will show a sharp drop at each transition point. It is expected that the polaron crossover will 
also show concomitant sharp jumps at each magnetic transition and the crossover to small 
polaron behavior will occur at lower value of g_ than in the FM limit. This is shown clearly in 
Fig. 4. It is further evident that in a double-exchange system, both the magnetic transitions 
and the electron-phonon coupling localize the electron, when the polaron crossover and 
magnetic transitions are overlapping. The locations of the large polaron region (A) and the 
small polaron region (B) are indicated in the vs J phase diagram (Fig. 1). In Fig. 1, 
for large values of (?_, the CAl and CA2 phases are very narrow and appear as a single 
phase boundary of FM-AFM region. So the line of separation of polaronic regimes (A and 
B) appears as a point for CAl and CA2 in Fig. I. Thus the results presented in Fig. 3 and 
Fig. 4 have a bearing on the transport behavior of our model. 

Having discussed transport we now redirect our attention to thermodynamic properties. 
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With this in mind we show in Fig. 5 the variation of the specific heat in the low temperature 
region in the FM state with zero and one phonon states. With application of an external 
magnetic field h, Cy takes lower values than for h = which is expected, as the average 
energy decreases with application of h in the FM state. For CA1(| Su \— 2), CA2(| Su |= 1) 
and AFM (| Su \= 0) states the external magnetic field will tend to align the core spins 
to feromagnetic order(| Su | = 3). For CAl, CA2 and AFM states at low field and low 
temperatures it can be shown from the present calculation that Cy does not change much 
from the h = limit as long as h does not shift | Su | to higher values. For larger h, 
as the ground state changes from lower | S12 | to a higher one, Cy decreases in the low 
temperature region. For CMR materials, there are some reports on measurements of Cy 



but these are measured in the FM state p2|, |13| , |14[ . In such cases, it was found that for low 
doping regions, Cy decreases with an increasing magnetic field. Our present calculation of 
Cy seemingly agrees with these experimental findings. The difference in the quantum and 
classical cases for specific heat, as far as the core spins are concerned, is exemplified in Fig. 
6 and Fig. 7 for FM and AFM cases respectively. The quantum results evidently yields the 
correct low-temperature limit. 



V. CONCLUSIONS 

The Anderson-Hasegawa model, though almost fifty years old, is able to capture all 
the crucial features of the double-exchange mechanism, originally proposed by Zener. The 
model is restricted to just two sites but the limitation should not be too serious when the 
electron hopping, influenced by thermal fluctuations, lattice distortions, phonon effects and 
other interactions, is expected to be incoherent. Incoherent hopping, albeit quantum in 
nature, is quite distinct from coherent band-like propagation, and approximately follows 
a Markovian process as far as the quantum diffusion of the electron is concerned. For a 
Markovian process only the pre-hopping and post-hopping sites matter. Therefore, the two 
site abstraction of the underlying three dimensional lattice provides the simplest paradigm 
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which can be exploited for analyzing a variety of phenomena which are of current interest 
in manganites. With this in mind, we have felt the need of carefully reexamining the exact 
quantum solution of the Anderson-Hasegawa model, for realistic values of the core spins. We 
have further used this model as the basic building block, in order to incorporate one or the 
other phenomena of relevance to manganites. These include superexchange and polarons, 
which have been the focus of our attention here. Indeed, we have found that Superexchange, 
when properly treated in conjunction with quantum spin dynamics of the core spins, leads 
to additional terms in the Anderson-Hasegawa Hamiltonian, which are absent in the clas- 
sical approximation of the core spins. A similar effect can also be observed with the exact 
solution of the Anderson-Hasegawa model, with the addition of phonon coupling. This may 
lead to the enhancement of the site energy term, which will have important consequences 
in these systems, especially in the paramagnetic state of the manganite systems. More- 
over, the discreteness associated with the effective hopping in this quantum case was shown 
to have further consequence for thermodynamic and transport properties. In conclusion, 
therefore, we find that the Anderson-Hasegawa model continues to remain relevant for the 
understanding of topically important issues in manganites. 
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Figure Captions : 

FIG. 1. The vs J phase diagram (h = 0) for | S ± | = | S 2 \= § and * = 1. (A) and (B) 
denote large polaron and small polaron region respectively. 

FIG. 2. Variations of ground state spin configuration | S\ + S 2 | with g_ for t — 1 and 
J = 0.01, 0.04 and 0.09, h = (in units of lo = 1). 

FIG. 3. Variations of t^ff and A corr / 1 g- with g^ for t = 1.0, J = 0.01 and /i = 0. 

FIG. 4. Variations of effective kinetic energy t%ff (dashed line) and polaron crossover 
^corr/9- (solid line) with g_, for t = 1.0, J = 0.09 and h — 0. The sharp jumps in and 
Korr/g- occur at values of where the magnetic transitions take place (see Fig. 2). 

FIG. 5. Variations of Cy (in arbitrary units) for g_ = 0.6, J = 0.01 and t — 1, for 
different values of the magnetic field /i = 0, 0.01, 0.05. 

FIG. 6. Variations of Cy (in arbitrary units) for g_ = 0.2, J = 0.02, h — and t = 1, in 
classical (solid line) and quantum (dashed line) formulation of the core spins. The ground 
state is FM. 

FIG. 7. Variations of Cy (in arbitrary units) for g_ = 0.9, J = 0.2, h — and t = 1, in 
classical (solid line) and quantum (dashed line) formulation of the core spins. The ground 
state is AFM. 
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